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Evanescent matter-waves produced by an atom wave packet incident onto a repulsive barrier edge 
can be back-refracted and reconstructed by the application of negative-index "comoving" potential 
pulses. One shows that those collapses and revivals of evanescent matter waves give rise to surface 
matter waves and should be observable via atom reflection echoes issued from the barrier interface. 
This property, together with the property of inducing negative refraction, makes such potentials the 
matter-wave counterpart of negative-index materials in light optics. 
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With the fast development of matter- wave optics, 
many of the functions previously operated in light optics 
have been realised: atom diffraction and mirrors, beam 
splitters, atom lasers, atom holography, etc. [1]. Specific 
characters of those processes originate in the properties 
of the associated particle: non-zero atom mass, vacuum 
dispersion for the "de Broglie" waves (implying longitudi- 
nal wave packet spreading) , scalar character of the atomic 
wave function, influence of the internal atomic degrees of 
freedom. Along this viewpoint, novel areas in the field of 
atom optics are presently explored, including e.g., the 
devising of non-diffracting atom nano-beams via a spe- 
cially designed transverse Stern-Gerlach interferometer 
[2] . Recently, we have proposed to extend the concept of 
"meta-optics" (known as negative index materials, NIM- 
s, in light optics [3-4]) to matter waves [5]. The main 
specificity of NIM "meta-media" for atom optics is their 
transient character, linked to the fact that, for matter 
waves, the reversal of atom group velocity is operated and 
can be only transient - contrary to the equivalent process 
in light optics where the phase velocity is reversed. In [5] , 
it was proposed to use a "comoving" potential pulse [6] 
to reverse the atom group velocity and induce negative 
refraction. In this letter, we examine the behaviour of the 
atomic wave packet in negative-index media and analyse 
its dynamics in classically forbidden regions (evanescent 
or semi-evanescent matter wave packets). The dynam- 
ics of matter waves impinging a static potential barrier 
is strongly altered by the comoving potential. In par- 
ticular, we investigate here the possibility of engineering 
evanescent matter-wave packets and predict the genera- 
tion of atom reflection echoes as well as surface matter 
waves at the barrier edge. 

It is commonly admitted that, in the semi-classical 
regime, light and matter waves involving a single direc- 
tion of space {e.g. x), behave identically provided they 
are monochromatic, i.e. provided they involve a single 
value of the wave number k (or particle's energy). For 
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both kinds of waves, the evanescent wave appears, at 
a normal incidence with respect to a planar interface, 
as soon as the refractive index of the second medium is 
imaginary. This means for light a conducting medium 
(metal, plasma). For particles, a constant potential en- 
ergy barrier Vq , extending from x = to x infinite, higher 
than the initial kinetic energy Eq of the particle, or, in 
the case of a wave packet, higher than the kinetic-energy 
distribution, provides an imaginary index since the index 
for matter waves in the region x > is n = (1 — Vo/i?o)^^^- 
In a 2D geometry (x, z) with the same potential barrier, 
the motion along z is free, which leads to a total reflection 
in the x < half plane and to an evanescent wave in the 
X > side. The corresponding index is static. However, 
as shown in [5] , a negative index for matter waves is nec- 
essarily a transient process, thus needing incident wave 
packets, i.e. some k- momentum distribution. As a con- 
sequence, light and matter evanescent waves behave dif- 
ferently as atomic wave packets do not factorize in x and 
t, due to vacuum dispersion. The existence of evanescent 
matter waves is one of the fundamental processes intro- 
duced by quantum mechanics, and has many outstanding 
consequences, mainly explored, up to now, for electrons 
(like electron tunnelling, as in tunnelling microscopy [7], 
or Josephson Effect in superconductivity [8] ) . For atomic 
systems, it has not been much studied, except in molec- 
ular spectroscopy [9], collision physics [10] and more re- 
cently cold atom physics [11], where evanescent waves 
are commonly involved {e.g. quasi-bound states and res- 
onances of various types). The evolution in x and t of 
a partially evanescent matter-wave packet, experiencing 
a pulse of comoving potential, will be described in the 
following. This comoving potential has the general form: 

V{t,x) ^ s{t) cos (2ttj^ (1) 

where s(t) is a signal restricted to a finite time interval 
{e.g. 0, Ti) and A is the spatial period. In the following 
we shall consider, as in previous studies [2,5], magnetic 
potentials acting on magnetic atoms, such as metastable 
argon atoms Ar*('^P2)- Actually any type of potential 
of the type (1) able to act upon the atomic wave packet 
motion would be convenient. Similarly, the repulsive po- 
tential barrier can be created using a static magnetic field 
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p(fc) being the momentum distribution and 




FIG. 1: Principle scheme of an experiment on evanescent and 
semi-evanescent wave packets submitted to comoving mag- 
netic potential pulses (see text). A static potential barrier is 
generated by the static magnetic field Bq. Potential pulses, 
comoving in the x direction, are produced by a generator ( "RF 
pulse"). The incident wave packet propagates in plane (x, z). 
The reflected wave packet and reflection echoes propagate in 
a direction of the (x, z) plane, symmetric of the incident di- 
rection with respect to the x axis. 



Bf) in the half-space a: > (see figure 1). 

We are dealing here with three wave functions, namely 
^i(t, x), the incident wave propagating from x = —oo to 
X = 0, ^r(i, x), the reflected wave propagating backwards 
from x = Otoa; = — oo, and the transmitted evanescent 
wave 5'(t, x) in the region a; > 0. For a specific value of 
the wave number k, i.e. a specific value E of the total 
incident energy, the three wave functions have the same 
time-dependence, in exp[—ihk'^t/ (2m)]. In the absence 
of comoving field, the conservation of the probability flux 
leads to the usual transmission and reflection factors (in 
amplitude): T{k) = 2k/ {k + in) and R{k) = T{k) - 1, 

{2m/h^)Vo > fc2. As ex- 
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pected, since the evanescent wave does not carry any 
probability flux, \B\ = 1. Similar expressions of R and 
T are obtained for k > a (partial propagation within the 
barrier), by replacing in by k' = Vfc^ — a^. In a first 
approximation, these expressions of R and T factors also 
hold if the comoving potential is applied at a small but 
finite distance from the barrier [0 < x << 1/k). In spite 
of the fact that, -since it does not propagate -, the evanes- 
cent wave never really leaves the barrier edge, we shall 
assume that the effect of the comoving potential pulse 
V(t, x) can be treated within the only half-space a; > 0. 
Under this condition, the time dependent Schrodinger 
equation is: 



ih a* 



2m 



a^* + [Vo + v{t, x)] ^ 



(2) 



where 'i/{t,x) is the wave function. Let '^Q{t,x) be the 
solution in absence of comoving potential, that is (up to 
a multiplicative constant): 
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Setting 



Mt,x,k) = e-^'^e-'^'''^''^'^ 



(j){t, x, k) — f{t, X, k)(l)o{t, X, k) 



(4) 



(5) 



and assuming that the "perturbation" factor f evolves, as 
a function of x, much slower than (/>o , then the term in 
d^f can be neglected with respect to ndxf ■ The validity 
of this latter approximation will be verified a posteriori. 
It leads for f to the equation: 



ihdtf - —ndj + Vf 
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Let us define two new (complex) variables 
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Using (7) it is readily verified that [dt + ^ 
With variables u, v, eq. (6) becomes: 

Setting / = A{u)exp[iip{u,v)], one gets d^ip 
and finally: 
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dv U(m, V ) 



(6) 
(7) 

) = 2d^. 

(8) 
-V/{2h) 

(9) 



A can be taken as a constant, in so far as its dependence 
on u is incorporated in the lower bound of the integral. 
For a comoving potential of the form (1), at a given value 
of u, setting v' — 2t' - u, one gets from (8): 
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where the upper bound is (t) — Min [t,Ti]. Finally, the 
phase shift depends on 3 independent variables, k (via 
k), t and x {via u). 

It is interesting to compare this expression of tp to that 
of the phase shift - let us call it (fpro - previously used in 
the case of a purely propagating matter wave (cf. Ref. 
[5], Eq. 8): 



{k,t) 



dt s(t) cos 



27r hk^ 
A m 



(11) 



Expression (10) becomes identical to (11) when in is re- 
placed by k, provided that |m| can be considered much 
smaller than \v\ . In other words, replacing in by k, this 
equivalence holds provided that |i — f|a:| << |i + ^a:|. 
This means that, in the case of propagation, all relevant 
values of t and x are supposed to be close to those related 
to the propagation of the wave packet centre. Such an 
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approximation is clearly not valid for an evanescent wave 
packet. The signal s(t) in Eq. (1) is taken as : 

s{t) = Ce^{t + e)"2e"*/^ /or < i < n ; = elsewhere 

(12) 

where C is a constant. For a magnetic potential C = 
g^sBmaxM , where g is the Lande factor, /is the Bohr 
magneton, B„iax the maximum value of the magnetic 
field magnitude and M the magnetic quantum num- 
ber. We consider here atoms the spin of which is J 
— 2, namely argon metastable atoms (Ar* ^P2) hav- 
ing a velocity along x axis of 4 m/s (de Broglie wave- 
length XdB = 2.8 nm), polarized in the M = +2 Zee- 
man sublevel. The maximum magnitude of the comov- 
ing field B^ax is varied from up to 400 Gauss. The 
centre of the Gaussian momentum distribution p{k) is 
ko = 2.244 lO^TO"^. Evanescent wave packets appear 
when k < ko within a repulsive potential barrier of mo- 
mentum height a = ko (this corresponds to a magnetic 
field of 127 Gauss). This value of the height has been 
chosen in order to get a "semi-propagative" regime of the 
wave packet. The standard deviation is Sk — 0.005 ko. In 
s(t), time parameters are different from those previously 
used since they must be adapted to typical characteristics 
of the wave packets in the region a: > 0. In the present 
case we take e = 7.4 /is, r — 0.37 /is, ri = 0.60 fis. The 
spatial period is A = 2 fim [12]. The effect of the co- 
moving pulse on the wave function in half space a; > is 
simply obtained by incorporating the factor exp{i(p) into 
the expansion of this wave function over k (Eq. 3). At 
this point, let us examine the validity of approximation 
used previously to get the simplified equation (6). Actu- 
ally it is rather poor for the original comoving potential 
in cos{2ttx/ A), in the sense that the ratio \k dxf \ / \d1f\, 
which is zero for fc = (fc = ko), becomes larger than 1 
only when \k — ko\ > 5k. The reason is that such a poten- 
tial starts abruptly at a; = 0. The situation is greatly im- 
proved by slightly shifting the spatial dependence, into, 
for instance, cos[2t:x/ K — tt/IO) which is quite easy from 
an experimental point of view. In such a case the ra- 
tio is larger than 10 as soon as |/c — /co| > Sk/2Q. In 
other words, only a narrow central slice of the spectrum 
violates the approximation. It is expected - and it has 
been verified - to be of little importance in so far as, for 
the spatial dependence of the wave packets, it involves 
large distances (x > 20/dk = 1.57 /im). In what follows, 
the spatial dependence in cos(2ttx/ A — tt/IO) has been 
adopted, at least when semi-evanescent wave packets are 
considered. For fully evanescent ones {i.e. a — 1.02 fco), 
this correction has not to be used in so far as the ques- 
tionable part of the spectrum (k « 0) corresponds to a 
very low intensity. 

The profile in t and x of the squared modulus jV'ol 
of the unperturbed evanescent wave is shown in figure 2. 
Figures 3 a-b show the squared modulus of the wave 
function modified by the comoving field, for B^ax — 50 
and 28 Gauss. Two effects are clearly seen in these fig- 
ures: first, a negative refraction effect bringing the wave 
packet backwards to the barrier edge after some delay; 




FIG. 2: Profile in t and x of the squared modulus \ipo\^ of the 
unperturbed semi-evanescent wave function, (a) The height 
a (in momentum) of the potential barrier is chosen equal to 
central value fco of the incident momentum distribution. At x 
— 0, the time dependence is Gaussian. For a; > 0, one observes 
as expected the decay characteristic of the evanescent part 
together with a partial wave packet propagating at a positive 
group velocity and broadening with time, (b) Same as (a) for 
a pure evanescent wave (a = 1.02 fco). No more propagation 
is observed: exponential spatial decay of the wave at t = 0; 
wave-packet fully vanishing for t > 0. 



simultaneously, a significant enhancement of the wave on 
its way back to the barrier edge. As the spatial range 
of the comoving potential is infinite, one could think the 
reflected part of the wave packet is also brought back 
to the barrier edge by the negative refraction. Actually 
this does not occur because the central momentum fco 
in the region a; < is much higher than any propagat- 
ing momenta within the barrier, the order of which is 
Sk. Consequently the magnitudes of Bmax considered 
here {Bmax — 50 Gauss) are not sufficient to bring the 
reflected part back on the barrier edge. Therefore the 
revival observed at a; = O'^ after some delay, originates 
from the matter wave inside the barrier. Because of the 
boundary continuity conditions at x = 0, this revival gen- 
erates a secondary wave packet, an atom "reflection echo" 
in the region a; < 0. Actually similar echoes are obtained 
even when the barrier height is such that the propagative 
part is largely suppressed, e.g. a = 1.02 fcp = fco + 4 5k. 
In such a case (see fig. 4a), no appreciable internal neg- 
ative refraction effect is seen {i.e., no propagation effect 
inside the potential barrier), whereas the wave echo re- 
mains. Note the narrowing of the echo (the width of 
which is similar to the initial one), produced by a time- 
reversal property of the meta-medium [13]. On the other 
hand, in the case of an almost pure propagation, e.g. for 
a = 0.99 fco) the negative refraction inside the barrier, 
accompanied by the narrowing of the wave packet [13], 
becomes the main phenomenon. Refiected and transmit- 
ted intensities emitted from the barrier are characterized 
by internal refiection and transmission factors given by 
expressions similar to the previous ones. The internal re- 
flection generates an additional wave packet which prop- 
agates towards positive x as soon as the comoving pulse 
is over. A similar enhancement effect is observed when a 
second pulse, delayed by 0.6 ^s, is appHed (see fig. 3c). 
This is a clear evidence for the beginning of a tempo- 
ral series of collapses and revivals of the semi-evanescent 
matter wave at the potential barrier. The same type of 
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FIG. 3: (Left side) Same as fig. 2a (semi-propagating wave 
packet, a = ko), in the presence of a single comoving potential 
pulse of duration ri = 500 ns and magnitude (see text) (a) 
Bmax = 50 Gauss, (b) Bmax = 28 Gauss; in (c), two identical 
subsequent pulses (Bmax = 28 Gauss) are applied, the second 
one being delayed by 600 ns with respect to the first one. Note 
that the vertical scale is magnified by 10. Sequential rebounds 
are seen. 

FIG. 4: (Right side) Same as fig. 3, with a barrier height 
a = 1.02 fco (pure evanescent wave). The comoving potential 
pulse has a magnitude Bmax ~ 35 Gauss, (a) Single comoving 
pulse of duration ri = 425 ns; (b) two sequential pulses are 
applied, within a time delay of 425 ns. This makes appear 
a second echo, (c) Evanescent wave revivals with a single 
potential pulse of duration 500 ns. 

result is obtained with a pure evanescent wave packet 
(see fig. 4b). Successive revivals of evanescent waves are 
also predicted when a single comoving potential pulse of 



sufficient duration is apphed (fig. 4c). This behaviour 
of the wave packet when subsequent pulses of comov- 
ing potential are apphed strongly suggests the existence 
of a surface matter wave, which remains confined along 
the barrier edge (x = axis) while exhibiting a series 
of temporal rebounds along the t axis. Because of the 
boundary conditions, each revival of the evanescent wave 
at the barrier generates in free space a new, retarded, 
refiected matter wave packet propagating in the a; < 
direction, an atom reflection echo which should be ob- 
servable. This should be a signature of the behaviour of 
the evanescent wave packet inside the barrier. Another 
observation mode would consist in cutting off the poten- 
tial barrier at some distance, smaller than the evanescent 
decay length, and measure the transmitted intensities. 

To summarize, adequate pulsed comoving potentials 
have proven to be able not solely to induce a negative re- 
fraction of matter waves, but also to cause drastic effects 
on the dynamics of evanescent or semi-evanescent matter 
waves, generating repeated enhancements in the vicinity 
of the barrier edge, which gives rise to surface matter 
wave, as well as atom reflection echoes. This makes these 
"negative index media" a counterpart of left-handed me- 
dia in light optics, with however important differences 
related to the different natures of matter and light waves. 
One should emphasize the close relationship between the 
present approach (i.e. interaction of virtual atoms of 
the evanescent matter- wave with "real" magnetic or light 
fields) and that demonstrated long time ago about the in- 
teraction of "virtual" photons with real atoms, like satu- 
rated absorption in evanescent light waves [14]. Various 
applications can be imagined for matter waves : some 
of them are similar to those of left-handed materials in 
light optics, such as atom wave focusing (meta lenses and 
atom nano- lithography) , cloaking, etc.; some others are 
specific of matter waves such as beam splitters and in- 
terferometers and guided matter wave along a potential 
barrier edge. 
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